HOLDER EQUICONTINUITY OF THE INTEGRATED 
DENSITY OF STATES AT WEAK DISORDER 
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Abstract. Holder continuity, \Nx{E) - Nx{E')\ < C\E ~ E'l"", 
with a constant C independent of the disorder strength A is proved 
for the integrated density of states N\ [E) associated to a discrete 
random operator H = Ho + XV consisting of a translation invari- 
ant hopping matrix Ho and i.i.d. single site potentials V with an 
absolutely continuous distribution, under a regularity assumption 
for the hopping term. 



1. Introduction 
Random operators on £^(Z^) of the general form 

= Ho + XV^ , (1) 

play a central role in the theory of disordered materials, where: 

(1) V^ip{x) = uj{x)ip{x) with uj{x), X G Z^, independent identi- 
cally distributed random variables whose common distribution 
is p{uj)duj with p a bounded function. The couphng A G M is 
called the disorder strength. 

(2) Ho is a bounded translation invariant operator, i.e., [S^, Ho] = 
for each translation S^iplx) = ip{x — ^ G Z*^. 

The density of states measure for an operator H^^ of the form eq. is 
the (unique) Borel measure dNx{E) on the real line defined by 

x:\x\<L 

and the integrated density of states Nx{E) is 

Nx{E) := / dNxie). 
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It is a well known consequence, e.g., ref of the translation invariance 
of the distribution of H^^ that the density of states exists and equals 

N^{E) = [ {60, P(-oo,i?)(i^J5o> dP(cu) , every E eR; 
Jn 

for P almost every uj, where P is the joint probability distribution for 
uj and Q = is the probability space. 

The density of states measure is an object of fundamental physical 
interest. For example, the free energy / per unit volume of a system 
of non-interacting identical Fermions, each governed by a Hamiltonian 
of the form eq. (|H), is 

= -Pj ln(l + e-^(^-'^))diV,(E) , 

where (3 is the inverse temperature and /i is the chemical potential. 
Certain other thermodynamic quantities (density, heat capacity, etc.) 
of the system can also be expressed in terms of A^a- 

Our main result is equicontinuity of the family {Nx{-), A > 0} within 
a class of Holder continuous functions, that is 

Nx{E + 6)- Nx{E -6) < Ca 5", for all A > , (2) 

under appropriate hypotheses on Ho- The exponent a < 1 depends on 
Ho as well as the probability density, with a = | at generic E for a 
large class of hopping terms if p is compactly supported. 

A bound of the form eq. for the integrated density of states 
associated to a continuum random Schrodinger operator is implicit in 
Theorem 1.1 of ref. although uniformity in A is not explicitly noted 
there. The tools of ref. PP carry over easily to the discrete context to 
give an alternative proof of eq. (0). However the methods employed 
herein are in fact quite different from those of ref. PP, and may be 
interesting in and of themselves. 

The main point of eq. (j2)) is the uniformity of the bound as A — 0, 
since the well known Wegner estimate fU], see also [7| Theorem 8.2], 

^Nx{E) IIpIL 

~Je~ - ^ ' 

implies that Nx{E) is in fact Lipschitz continuous, 

Nx{E + 6) - Nx{E - 6) < M- 25 . (4) 

However, the Lipschitz constant \\p\\^ /A in eq. (jH) diverges as A ^ 0. 
Such a singularity is inevitable for a bound which makes no reference to 
the hopping term, since dNx{E) = X~^p{E/ \)dE for Ho = 0, as may 
easily be verified. However if the background itself has an absolutely 
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continuous density of states, the Wegner estimate is far from optimal 
at weak disorder. 

The translation invariant operator Ho may be written as a superpo- 
sition of translations, 



where 



dq 



is the inverse Fourier transform of a bounded real function e on the 
torus T"^ = [0, ^nY, called the symbol of Ho- For any bounded mea- 
surable function /, 

dq 



{2-kY_ 



from which it follows that the density of states No{E) for Ho obeys 

dq 



J f{E)dNo{E) = J^J{e{^)) 



(27r)' 



In particular. 



We define a regular point for £ to be a point G M at which 

No{E + 6) - No{E - 6) < T{E)6, (5) 

for some T{E) < oo. In particular if e is and Ve is non-zero on the 
level set {^(q) = E}, then is a regular point. For example, with Ho 
the discrete Laplacian on 

Ho^pix) = ^{x + I) + tp{x - I) , 

we have the symbol 6{q) = 2cos(g) and every E G (—2,2) is a regular 
point. However at the band edges, E = ±2, the difference on the right 
hand side of eq. (0) is only 0(62), and these points are not regular 
points. We consider the behavior of Nx{E) at such "points of order 
a," here a = 1/2, in Theorem 01 below. 

Our main result involves the density of states of Hx at a regular 
point: 

Theorem 1. Suppose J \uj\'^ p{uj)duj < 00 for some 2 < q < 00 or that 
p is compactly supported, in which case set q = 00. If E is a regular 
point for e, then there is Cq = Cq{p, T{E)) < 00 such that 

Nx{E + 6) - Nx{E ~ 6) < r{E)6 + A^^^+i^(5^(^~t^ (6) 
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for all X,S>0. 

For very small S, namely 

A ~ 

the Wegner bound eq. Q is stronger than eq. Thus Theorem [T] 

is useful only for 

29+1 
5 > Xi+^ . 

Combining the Wegner estimate and Theorem^ for these separate re- 
gions yields the following: 

Corollary 2. Under the hypotheses of TheoremUi there is Cq < 00, 
with Cq = Cq{p, T{E)), such that 

Nx{E + 6) ~ Nx{E - 6) < Cq6^^^'^^ (7) 
for all X,6>0. 

Thus, the integrated density of states is Holder equi-continuous of order 
i as A — >• (if p is compactly supported). 

The starting point for our analysis of the density of states is a well 
known formula relating dA^A to the resolvent of H^^, 



dE 



lim - [im {So, (H^-E- ir/)-'(5o> 



00) 



The general idea of the proof is to express Im ((5o, [H^ — E — \r]) ^5o) 
using a finite resolvent expansion to second order 

{H^-E-\ri)-^ 

= {Ho-E- ir])'' - X{Ho ~E- ir])-W^{H, -E- it])-' 
+ X\H, -E- ivr'V^iH^ -E- iv)-'V^{Ho -E~ ir^)-^ , 

and to use the Wegner bound eq. (jH)) to estimate the last term, with 
the resulting factor of 1/A controlled by the factor A^. 

Here is a simplified version of the argument which works if E falls 
outside the spectrum of Ho and ipE = {Ho — E)~^6o G £^{Z'^). The first 
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two terms of eq. (jHI) are bounded and self-adjoint when rj = 0, so 

lim - [ Im (So, {H^-E- i7]Y^5o) dP(cu) 

= A^lim- [ lm{^E,V^{H^-E-ir])-'V^^E)dnuj) 
'7io Jn 



< 



xHimJ2\Mx)\\My)\ 



\ujixMy)(^5,,{iH^-Ey + r]y'6y)\dnuj). 
If p is, say, compactly supported, then 



lim — 

J?iO TT 



u;{x)u{y) (5^, {{H^ - Ef + rf) 5y) dP(^) 

lim^ / l5,,{{H^-Ef + rfYUy)d¥{Lo) 



< 



< i 



by the Wegner bound, and therefore 
dN^{E) 



dE 



< A 



||2 

eWi , 



for E ^ a{Ho) 



(9) 



We have used second order perturbation theory to "boot-strap" the 
Wegner estimate and obtain an estimate of lower order in A. Un- 
fortunately, as p was assumed compactly supported, E is not in the 
spectrum of Hx for sufficiently small A, and thus dNx{E) / dE = 0. So, 
in practice, eq. © is not a useful bound. 

Nonetheless, in the cases covered by Theorem ^ Hx can have spec- 
trum in a neighborhood of E, even for small A, since E may be in 
the interior of the spectrum of Ho- Although, the above argument 
does not go through, we shall exploit the translation invariance of the 
distribution of H^^ by introducing a Fourier transform on the Hilbert 
space of "random wave functions," complex valued functions \E'(a;, cj) 
of (x,tu) e £2(Z^) X n with 



X 



UJ) < oo . 



Under this Fourier transform an integral of a matrix element of 
f{Hij) is replaced by an integral /^^ over the (i-torus of a matrix element 

of f{H]f), with Hk a certain family of operators on (see eq. (UHl)). 

Off the set Se ■= {k G T'^||£:(k) - E\ > e} with e >> 6, we are able to 
carry out an argument similar to that which led to eq. Q. To prove 
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Theorem ^ we shall directly estimate 

N{E + 6)-N{E-6) = [ {6o,Ps{H^)So)dF{uj) , 

Jn 

with Ps the characteristic function of the interval [E — 6,E + 6], because 
the integrand on the r.h.s. is bounded by 1. Since E' is a regular point, 
the error in restricting to S'^ will be bounded by r{E)e. Choosing e 
optimally will lead to Theorem ^ 

More generally, we say that ii^ is a point of order a for e, if there 
exists T{E] a) such that 

No{E + 6) - No{E - 6) < r(E;a)r. 

U E ^ a{Ho), we say that E is a. point of order oo and set T{E; oo) = 0. 
For points of order a we have the following extension of Theorem ^ 

Theorem 3. Suppose J \uj\'^ p{uj)duj < oo for some 2 < q < oo or that 
p is compactly supported, in which case set q = oo. If E is a point of 
order a < oo for e, then there is Cg^a = Cg^a{p,^iE; a)) < oo such 
that 

NxiE + 5) - NxiE - 5) < T{E;a)5" + 
for all A,5 > 0. 

When a = oo and q = oo, so E ^ cr^Ho) and p is compactly sup- 
ported, the result is technically true but uninteresting since E ^ cr{H\) 
for small A, as discussed above. However for g < oo, we need not have 
that p is compactly supported, and E ^ cr{Ho) may still be in the 
spectrum of H\ for arbitrarily small A. In this case, eq. (fTUI) signifies 
that 

NxiE + 5) - NxiE ~ 5) < Cg,^X'^U'~l , 

which in fact improves on the Wegner bound for appropriate A, 6. 

As above, we may use the Wegner bound for 6 very small to improve 
on eq. (fm|l : 

Corollary 4. Under the hypotheses of Theorem \^ there is Cg^a = 
Cq,a{p,^{E] a)) < oo such that 



Nx{E + 6) - N^{E - 6) < Cg,J 
for all A,5 > 0. 

The inspiration for these results is the (non-rigorous) renormalized 
perturbation theory for dA^A which has appeared in the physics lit- 
erature, e.g., ref. |Hj and references therein. If J ujp{ijj)du = and 
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J u'^p{uj)duj = 1, as can always be achieved by shifting the origin of 
energy and re-scahng A, then the central result of that analysis is that 



dE TT \ 

where Tx{E), the so-called "self energy," satisfies lmTx{E) > with 

limlmrA(^) ^ limIm(5o, (//o-^-ir/)"^5o) = tt^^^^ . 

Up to a point, the self-energy analysis may be followed rigorously. 
Specifically, one can show (see 

Proposition 1.1. // J Lop{uS)duj = and J uj'^p{uj)duj = 1, then for 
each A > there is a map Fx from {Imz > 0} to the translation 
invariant operators with non-negative imaginary part on £^(Z^) such 
that 

[ (H^-zy'dFico) = {H,-z-\^Txiz))-' , (11) 
Jn 

and for fixed z G {Im^; > 0} 

lim (5,, Tx{z)6y) = {6o, {Ho - z)-'6o) 6^,y . (12) 

However there is a priori no uniformity in z for the convergence in 
eq. (fT^ . so for fixed A we may conclude nothing about 



lim (Ho-E-m- \'^Tx(E + in)) 



-1 



Still, one is left feeling that Theorem ^ and Corollary |^ are not- 
optimal, and the "standard wisdom" is that something like the follow- 
ing is true. 

Conjecture 5. Let p have moments of all orders, i.e., J \uj\'^p{uj) < oo 
for all q > 1. Given Eo G M, if there is 6 > such that on the set 
{q : |£:(q) — Eo\ < 6} the symbol e is with V£(q) ^ 0, then there is 
Cs < oo such that 

for allXeR and E e [E^ - ^6, Eo + |5] . 

Remark: The requirement that p have moments of all orders is simply 
the minimal requirement for the infinite perturbation series for {Ho — 
z — AV^)~^ to have finite expectation at each order (for Im^; > 0). In 
fact, this may be superfluous, as suggested by the example of Cauchy 



JEFFREY H SCHENKER 



randomness, for which the density of states can be exphcitly computed, 
see ref. [7]: 

^Nx{E) = - [ i,, , , for piuj) 



71 Jt^ {e{q) - EY + A2 {2'kY ' '^^ ' tt 1 + cj2 
although J p{uj) \uj\'^ = oo for every q > 1. 

2. Translation invariance, augmented space, and a 
Fourier transform 

The joint probabihty measure F{uj) for the random function u : Z"^ 



IS 



fd 



on the probabihty space VL = W' . Clearly, P( uj) is invariant under the 
translations : VL ^ VL defined by 

^uj{x) = uj{x — ■ 

In particular, since 

S^Hi^jS^^ = Ho + Vri^uj = Hti-U) , (13) 

and S^H^S^^ are identically distributed for any ^ G Z"^, 
To express this invariance in operator theoretic terms, we introduce 
the fibred action of if^ on the Hilbert space L'^{Q; ^^(Z'^)) - the space 
of "random wave functions" - namely. 

We identify ^^(fi; f{Z'^)) with L'^{n x Z^) and denote the action of 
on the latter space by H, so 

[H^](a;,x) = ^e{^)'^{u,x - ^) + Aa;(x)^(a;, x) . 

The following elementary identity relates f{Hi_j)dF{uj) to /(H), for 
any bounded measurable function /, 

/ dF{u){6,J{H^)6y) = (Et4,/(H)Et5,) , (14) 
Jn 

where E^" is the adjoint of the linear expectation map K : L'^{Qx 7/) — > 
£2(Z'^) defined by 



[E^](x) = / m{uj,x)d¥{uj) 



Note that E"'' is an isometry from ^^(Z'^) onto the subspace of functions 
independent of u; - "non-random functions." 
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The general fact that averages of certain quantities depending on if^ 
can be represented as matrix elements of H is known, and is sometimes 
called the "augmented space representation" (e.g., ref. jHElEl) where 
"augmented space" refers to the Hilbert space L^(f2 x There are 
"augmented space" formulae other than eq. (HD), such as 

/ &^{uj)uj{x)uj{y)l,b.,,i{Ej)by) = (Et4,V/(H)VEt5^) , (15) 
Jn 

and 

/ d¥{u){S,,f{H^)So){So,g{H^)Sy) = (E^^, /(H)Po^7(H)Et5,> , 
Jn 

where Pq denotes the projection Pq'^^UjX) = ^'(a;,0) if x = and 
otherwise. The first of these (eq. p5j) ) will play a roll in the proof of 
Theorem Q 

There are two natural groups of unitary translations on L'^{Q x Z'^): 
and 

Note that these groups commute: [S'^,Tg'] = for every G Z'^. A 
key observation is that the distributional invariance of H^, eq. (fT^. 
results in the invariance of H under the combined translations T^S^ = 

S^T^HTIsI = H. 

In fact, let us define 

Ho = V^{uj,x) = uj{x)^{uj,x) . 

Then H = Ho + AV where Ho commutes with S'j and while for V 
we have 

To exploit this translation invariance of H, we define a Fourier trans- 
form which diagonalizes the translations S^T^ (and therefore partially 
diagonalizes H). The result is a unitary map : L^(f2 x Z*^) 
L'^{Q X T'^), with T'^ the c?-torus [0,27?)*^. Let us define J" first on 
functions having finite support in Z'^ by 
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It is easy to verify, using well known properties of the usual Fourier 
series mapping ^^(Z'^) — > L^(T''^), that T extends to a unitary isomor- 
phism L^{^ X if) ^ L2(Q X T^), i.e. that JTjTt = 1 and T^T = 1 
where T'^ is the adjoint map 

Another way of looking at T is to define for each k e T*^ an operator 
JTk : L2(n X Z'^) ^ L2(0) by 

JTk* = lim y e"^*^-^ J^.T.* , 

L— >oo ' 

where J is the evaluation map J^{ijj) — ^{oj, 0). The maps are 
noi bounded, but are densely defined with J^k^ £ L^(Q) for almost 
every k, and 

:^^'(cj,k) = J^k^(^) a.e. u;, k. 
If we look at L'^{VL x T*^) as the direct integral J® dkL2(n), then 

® 

dkJ'^k . 

This Fourier transform diagonahzes the combined translation S^T^^, 

as follows from the following identities for S and T, 

where, on the right hand sides, denotes the operator T^ipiuj) — 
iP{t^^u!) on L^(0). Furthermore, explicit computation shows that 

J^kV = u;(0)JP-k , 

where ci;(0) denotes the operator of multiplication by the random vari- 
able a;(0), ^ u;{0)'ip{u;). Putting this all together yields 

Proposition 2.1. Under the natural identification of L'^{Q,T'^) with 
the direct integral J dkL^(Q), the operator H = JTIJ^^ is partially 
diagonalized, H = Hyg^, with ifk operators on L^(Q) given by the 
following formula 

= ^e-*-«£(-e)T^ + Ac.(0). 




DENSITY OF STATES AT WEAK DISORDER 



11 



c c UT'^ 



(2vr)^ 



Let us introduce for each k G T'', 

dq 

so Hk = + Xuj{0). Note that 

where = 1 for every G That is, is an eigenvector for 

Applying the Fourier transform to the right hand side of the "aug- 
mented space" formula eq. ()14p we obtain the following beautiful iden- 
tity, central to this work: 

dk 

Similarly, we obtain 



u){5,JiH^)5y) = J^^^^e^<^-y)^x^,f(H^)x^) . (16) 



d¥{uj)uj{x)uj{y) {6^J{H^)6y) 

dk 



(277)'^^'^'^'''^ V^^^^^^ /(//k)^(O)xn) (17) 

from eq. (|15|). Related formulae have been used, for example, to derive 
the Aubry duality between strong and weak disorder for the almost 
Mathieu equation, see ref. |2] and references therein. 

As a first application of eq. (fT^ . let us prove the existence of the self 
energy (Prop. II. Ij) starting from the identity 

Proof of Prop. \1.1[ Since xn is an eigenvector of and 

{Xn,uj{0)Xn) = j t^p(c^)da; = 0, 
the Feschbach mapping implies 

(xu\h^-zY\^ = {e{k)-z-\^V^{z-'k)Y\ (18) 

^In fact, if £ is almost everywhere non-constant (so Ho has no eigenvalues) then 
e{k) is the unique eigenvalue for and the remaining spectrum of is infinitely 
degenerate absolutely continuous spectrum. One way to see this is to let (/>„(«) be 
the orthonormal polynomials with respect the weight p(f), and look at the action 
of on the basis for L?{Vl) consisting of products of the form YixeZ'' 't'n{x){^{^)) 
with only finitely many n(x) ^ 0. 
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with 

where P"*" denotes the projection onto the orthogonal complement of 
Xn in L\n). 

Let the self energy T\{z) be the translation invariant operator with 
symbol Tx{z; k), i.e., 

{S.,T.(^)S.) = X/"-'-"'r.(.;k)p-^. 

Clearly Tx{z) is bounded with non- negative imaginary part. Further- 
more by eq. (fTHj) and eq. (fTS|) . the identity eq. (fTT|) holds, namely 



[ {H^ - zyhlFicu) = {H,-z-X^T,iz)) 
Jn 



-1 



It is clear that 

from which eq. ()12|1 follows easily. □ 

3. Proofs 

We first prove Theorem ^ and then describe modifications of the 
proof which imply Theorem El 

3.1. Proof of Theorem 12 Fix a regular point E for e, and for each 
(5 > let 



fs{t) = {X{E-S,E+S)i't) + X[E-5,E+5]{t)) 



1 , t,e{E-5,E + 5) 

\ t = E±5, 

t^[E-5,E + 5]. 



Since N\{E) is continuous (see eq. (jH)), 

Nx{E + 6) - Nx{E - 6) = [ {6oj5{H^)So) dF{uj) 

Jn 

Thus, in light of eq. (fTB|) . our task is to show that 



'J'd 



dk 



Xn, fs{H^)Xn) 7^ < r{E)6 + A^^^+l^^^^^^ , (19) 



(2vr 
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with a constant Cq independent of 5 and A. Note that for each k e T'^ 



< 1 



^(i+l)Aia-f) 



so we can afford to neglect a set of Lebesgue measure A^*' 
on the l.h.s. of eq. (fT^. 

Consider k G T'^ with |e(k) -E\> 5. Then 



h{K)Xn = fsieik))xn = . 



Thus 

Xn, f5{H]i)xn 



E+S 



hm - / lra{ Xn 

V^O IT Je-5 



Xn,{fs{Hu)-fs{H^))Xn 
1 1 



A hm — 



E+S 



Im- 



Hi^-t-ii] H^-t -17] 
1 / 1 



Xn ) dt 



E-& t + ir]- eik) \^^' H^-t-irj 
1 



^(0)xq ) dt 



(20) 



since (t — £:(k)) ^ is continuous for t G [E — 5, E + 5]. Using again that 
fs{H^)Xn = 0, we find that the final term of eq. piU) equals 



- 5(k) 

E+S 



- e(k; 



■AM 



'^^O TT t - e(k) t + 17] - £(k) 



X 



(21) 



A(^(0)xn,-^^^t^^(0)xn 
Putting eqs. ((20)) and (PT|) together yields 



Xn,A(^^k)xn) = A^(^(0)xn, ./'™... o;(0)xn 



(^k-e(k))^ 



< A^ 



(Kk)-i5;|-5)^ 



^(0)xn,/5(^k)^(0)xn 
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Thus, for any e > 5, 



{\e(k)-E\>e} 



(e - 5)2 

where in the last equahty we have inverted the Fourier transform, using 
eq. ()17p . We may estimate the right hand side with Holder's inequality 
and the Wegner estimate: 



1- 



< ||cu(0)|P(^25 



A 



1-^ 

<3 



since {Sq, fs{H^)Sof < {So, f5{H^)So) for p > 1 (because {Sq, fs{H^)So) 
< 1). Here ||cj(0)||^ = /cu(0)«dP(w) for g < oo and ||cj(0)||^ = 
ess-sup^ 11^(0) I- 
Therefore 



Xn,fs{H^)Xn) < nE)e + A^^^^ 11^(0) ||J 25 

(22) 

where the first term on the right hand side is an upper bound for 

dk r dk 



2 



1- 

9 



/* dk /* 

/ (xn, f5iHk)xn) jTT^ < I 

J \\e(V.)-E\<e} \ / [^^) J\k 



'{|e(k)-£;|<e} ^ ' K^l^j J {\eO<.)~E\<e} 

Upon optimizing over e G (5, oo), this implies 
Jn 

which completes the proof of Theorem ^ D 

3.2. Proof of Theorem [31 If instead of being a regular point, E' is a 
point of order a then the proof goes through up to eq. , in place of 
which we have 



{xn, fsiH^)Xn) < r(E; «)e°+A2— ^ ||a;(0) g 



1- 
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Setting 6 = 6 + X'^'d' and choosing 7, l3 sucli tliat the two terms are of 
the same order yields 

2 + a \ qj 2 + a \ qj 

which imphes 

Jn 

completing the proof. □ 

Acknowledgements. I thank G.M. Graf and M. Disertori for stimu- 
lating discussions related to this work. 
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